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Abstract 

The multiplication rates of hot carriers in CdSe quantum dots are quantified using an atom- 
istic pseudopotential approach and first order perturbation theory. Both excited holes and elec- 
trons are considered, and electron-hole Coulomb interactions are accounted for. We find that 
holes have much higher multiplication rates than electrons with the same excess energy due to 
the larger density of final states (positive trions). When electron-hole pairs are generated by 
photon absorption, however, the net carrier multiplication rate is dominated by photogenerated 
electrons, because they have on average much higher excess energy. We also find, contrary to 
earlier studies, that the effective Coulomb coupling governing carrier multiplication is energy 
dependent. We show that smaller dots result in a decrease in the carrier multiplication rate 
for a given absolute photon energy. However, if the photon energy is scaled by the volume 
dependent optical gap, then smaller dots exhibit an enhancement in carrier multiplication for a 
given relative energy. 
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Solar cells normally produce a single electron-hole pair per photon absorbed. The excess pho- 
ton energy hco — E gap (where (0 is the photon frequency and E gap is the semiconductor band gap) 
is converted to heat, and represents a net loss for the photo-conversion efficiency. It is theoretically 
possible, however, for a photon of sufficiently high energy to generate two or more lower energy 
excitons, in a process known as multiple exciton generation (MEG). Provided that this process 
occurs faster than competing carrier relaxation processes, more than one pair of charge carriers 
might be collected per photon absorbed. This raises the prospect of designing solar-cell devices 
that exploit MEG in order to utilize high energy photons with greater efficiency than is currently 
possible. 

Early work in this area^ suggested that the MEG efficiency is enhanced by quantum con- 
finement of carriers, motivating experimental investigations of MEG in different nanostructured 
materials. Initial reports of very high MEG yields in PbSe, PbS, and CdSe nanocr 
were followed by a critical analysis of the assumptions used to extract the MEG yield from optical 
spectroscopy measurements ,L1511 6|2 1 1 24 | conse q uences f the energy dependence of the absorp- 
tion cross section, lS the non-linear dependence of the bleaching ratio on the number of excitons,^ 
and the non-radiative recombination of charged excitons^ on the accuracy of MEG measurements 
have all been discussed in the literature. This has led to a re-assessment of earlier experimental 
results. 

A controversy has recently emerged, though, as to how much, if at all, quantum confinement 
really enhances MEG. Most experimental investigations of MEG in nanostructures have focused 
on the MEG quantum yield (QY), defined as the average number of excitons produced per absorbed 
photon. MEG manifests itself by producing more than one exciton (QY>1) when the photon energy 
exceeds a certain threshold, hco t h. Experiments for PbSe 3 and Si^ nanocrystals have revealed 
that the QY measured at a given absolute photon energy decreases when the nanocrystal size 
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decreases, while other experiments have reported a weak and/or non-monotonic dependence of 
the QY on size.^1 When measured with respect to the relative photon energy, hco/E gap , however, 
the QY increases as the nanocrystal size decreases.^ 10 ' 1 1 1 16 1 Rabani et al. 1 ^ calculated the QY of 
a set of small CdSe and InAs nanocrystals and found an increase as the size decreases. A few 
experimenta l 1 ° l 14 l S I and theoretical^ 19 * 26 ! investigations have also attempted a comparison of the 
QY of quantum dots vs. bulk semiconductors. Beard et al.^reported a significantly higher QY in 
Si nanocrystals than in bulk Si at the same relative photon energy. Nair et al.^ found that the QY 
of PbS nanocrystals is lower than that of PbS cystalline films 27 when measured at fixed absolute 
photon energy. Similar results where reported by Pijpers et al.J^Delerue et al.^ and McGuire 
et alJ22 for PbS and PbSe nanocrystals. Notwithstanding these results, systematic experimental 
or theoretical studies of the size dependence of MEG are not available. Furthermore, the QY 
depends on to the ratio of the MEG rate, Tmeg, an d the total exciton decay rate, T tot . The latter 
includes competing decay processes, such as phonon-assisted carrier cooling. Since both Tmeg 
and T tot depend on the nanostructure size as well as the photon energy, the observation of a certain 
dependence of the QY on size does not necessarily imply that the same trend exists for Tmeg- 
Thus, it would be very useful to quantify the size dependence of the MEG rate itself. This quantity 
is difficult to measure experimentally, but accurate calculations of Tmeg as a function of size can 
provide information on multi-exciton dynamics and clarify the role of quantum confinement in the 
MEG process. 

Results and Discussion 

We use atomistic pseudopotential calculations, coupled with first-order perturbation theory, to 
quantify the dependence of the MEG rate on photon energy and nanoparticle volume for zinc- 
blende CdSe nanocrystal quantum dots (QDs). The volume dependence of the MEG rate is de- 
composed to show the relative contributions of Coulomb matrix elements and final density of 
states. Both hole and electron contributions to the MEG rate of an electron-hole pair are accounted 
for, so that we can make clear which of these dominates MEG. We find that while the hole MEG 
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rate is significantly larger than the electron MEG rate for a given excess energy, the electron decay 
dominates the MEG of photo-generated electron-hole pairs because of the larger excess energy of 
the electron. Tables 1-3 describe the energy and size dependence of the electron, hole, and exciton 
MEG rates of CdSe QDs. We find that (i) the MEG rate increases rapidly with the excess energy of 
the carriers, (ii) At fixed absolute photon energy, the MEG rate of a photo-generated electron-hole 
pair decreases as the volume of the QD decreases. However (iii) at fixed relative photon energy, 
the MEG rate increases as the volume decreases. 
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Figure 1 : Schematic diagram of MEG decay paths for an excited electron (left) and an excited hole 
(right). In practice, all possible combinations of electrons and holes in the final trion states are 
considered in the calculation of MEG rate for a given initial, excited state. 

When a carrier has sufficient excess energy (larger than ~ E gap ) it can spontaneously decay 
into a trion in a process known as impact ionization (see Fig. 1). The impact ionization rates for 
electrons and holes are calculated here by applying Fermi's golden rule to the formation of trion 
states. We consider here four nearly spherical QDs: Cd2i6Se2i3 (diameter D=2.9 nm), Cd3i2Se32i 
(D=3.2 nm), Cd484Se495 (D=3.7 nm), and Cd784Se739 (D=4.3 nm). Figure 2 shows the calculated 
MEG rates of electrons and holes for the two CdSe QDs considered here (diameter 2.9 and 3.7 nm, 
respectively). The MEG rates are plotted in Fig. 2 as a function of the absolute electron/hole energy 
measured with respect to the vacuum level (Figs. 2a and 2b), as well as the relative electron/hole 
excess energy \E — EQ\/E gap (Figs. 2c and 2d). Here Eq is the conduction-band minimum (CBM) in 
the case of electrons and the valence-band maximum (VBM) in the case of holes, and the optical 
gap Eg ap is calculated as E gap = £cbm — £vbm — Jvbm,cbm- Fig. 2 shows that the MEG rate of 
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electron and holes increases rapidly as their excess energy increases. Fig. 2 also shows that, for the 



same relative excess energy, excited holes have transition rates more than an order of magnitude 



greater than excited electrons. 




Figure 2: Calculated MEG rate as a function of the electron and hole energy for CdSe QDs with 
diameters of 2.9 nm (left) and 3.7 nm (right). Energies are shown in two ways, (a) and (b): 
with respect to the vacuum level; (c) and (d): with respect to the band-edge energy Eq (CBM 
for electrons and VBM for holes) and rescaled by the optical band gap E gap . The small symbols 
correspond to the transition rate of individual electron (blue) or hole (red) initial states to final 
trion states, while the large symbols in the bottom figures represent the arithmetic average of initial 
states within a 100 meV energy window. Solid lines, T(E), are fitted using the averaged data (large 
symbols). Arrows indicate the positions of the CBM and VBM. 



The calculated MEG rates for all the QD sizes considered here are summarized in Fig. 3a. For 
each QD, the energy dependence of the electron and hole MEG rates were fitted to a power law, 
Th e {E) oc [E — E t h) a ' ie , where the threshold energy E t h is the ground-state energy of a negative 



trion (for electrons) or a positive trion (for holes), calculated according to [equation] [8] [] 8 and 



[equation] [9] []9[ respectively. Table 1 summarizes the values of CC/, and Ct e extracted from the fit. 

We find that a varies between 2.1 and 2.7, depending on the QD size, and is slightly larger for 

holes than for electrons. The volume dependence of the MEG rates was calculated for a set of 

given electron and hole energies and fitted to the function Yj i e {V) = rjj g +Bf l>e V^ k < e , where e , 

e and % if , are fitting parameters listed in Table 2. In particular, /3/, e ranges between -2.0 and 

-2.4 for holes (depending on their energy) and between -1.2 and -1.8 for electrons. 

Table 1: Fitted parameters for the energy dependence of electron and hole MEG rates: T^ e (E) = 
Ahe(E — E t h) ah > e as shown in Figure 3a. 



Diameter (nm) 




oc e 


2.9 


2.3 


2.1 


3.2 


2.4 


2.3 


3.7 


2.7 


2.5 


4.3 




2.6 



Table 2: Fitted parameters for the volume dependence of electron and hole MEG rates: F^ e (V) = 
e + Bh,eV^ h e at given carrier energy, E. 



E(eV) p h B h 


E(eV) p e B e 


-7.8 -2.0 -11.2 
-8.0 -2.4 -66.8 
-8.2 -2.2 -69.7 
-8.4 -2.3 -128.0 


-0.6 -1.8 -0.93 
-0.7 -1.7 -0.78 
-0.8 -1.6 -0.60 
-0.9 -1.2 -0.18 



In order to explain the trends displayed in Fig. 3a, we decompose the MEG rate into the trion 
density of states (DOS) contribution, p(E) = Y,f8(E — Ef), and the effective Coulomb coupling 
contribution, \W\ 2 = KT(E)/2np(E). These two components are shown in Figs. 3b and 3c, re- 
spectively. The trion DOS (Fig. 3b) shows an energy dependence very similar to the overall MEG 
rate of Fig 3a, i.e., the number of final trion states increases rapidly as the excess energy increases. 
While this trend holds for trions deriving from both holes and electrons, the positive trion DOS is 
significantly larger than the negative trion DOS, as shown in the inset plot of Fig. 3b for a 3.7 nm 
dot. The reason for this behavior is that the manifold of hole states is much denser than its electron 
counterpart. As a result, the number of available trion states, for a given excess energy, is much 
higher for holes than electrons, and this dominates the overall MEG rates. Previous calculations 
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using an atomic pseudopotential approach for PbSe QDs have revealed a dense feature of valence 
states. 35 In the case of CdSe QD, the asymmetry in the density of states between the conduction 
and valence states is even more pronounced. 

The effective Coulomb coupling, \W\ 2 , (Fig. 3c) peaks just above the threshold energy and then 
exhibits a decay as the energy increases. The decay becomes faster as the size of the QD increases. 
The size dependence of the effective Coulomb coupling was obtained by fitting the data for all 
four CdSe dots to a power law: \W\ 2 °c V y . The Coulomb coupling scaling is energy dependent 
with / ranging from -2.6 to -3.2 for electrons and -3.6 to -4.0 for holes. This trend holds for all 
four QDs and is in contrast to the assumptions and conclusions of earlier studies. Specifically, a 
constant Coulomb coupling was assumed for an examination of impact ionization in bulk Si,^ 
and a recent tight-binding study concluded that there is only a weak energy dependence in the 
Coulomb coupling for PbSe QDs. 36 Recent pseudopotential calculations for CdSe QDs 25 found 
that the Coulomb coupling has no energy dependence for dots up to 3.0 nm. Our calculations 
suggest that this is not the case for larger QDs. In any case, the downward trend in the energy 
dependence of the Coulomb coupling is overshadowed by the opposing trend in the trion DOS. 

To understand the decrease in the effective Coulomb coupling with increasing excess energy, 
we analyze the spatial distribution of the initial-state electron and hole wave functions. For each 
single -particle state i of a QD, we calculate the fraction of the wave function norm that resides in 
the surface region of the QD, 5/ = J reS | V/(r) \ 2 dr, where 5 is a thin shell around the surface of the 
QD. Figure 4 shows S,- (Fig. 4a) and the charge density distribution of hole states located between 
-8.0 eV and -7.5 eV (Fig. 4b) for the 3.7 nm CdSe QD. The plot of 5 ; (Fig. 4a) indicates that states 
near the band edges are localized in the core of the QD. As the energy level moves away from 
the band edges, though, the wave functions spread to occupy a large portion of the surface shell. 
The broad peak located 3 eV below the VBM originates from surface states associated with the 
ligands used to passivate the QD surface (Fig. 4b). Also shown in Fig. 4a is the Coulomb coupling 
\W\ 2 of the same 3.7 nm QD (already shown in Fig. 2b). In the case of the electron levels, \W\ 
decreases as the surface projection, Su increases. However, this is not the case for the hole levels 
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Figure 3: (a) MEG rate T(E), (b) trion density of states p tr i 0n (E), and (c) effective Coulomb 
coupling, \W\ (E), as a function of the electron and hole energy for four CdSe QDs. Insets in (b) 
and (c) show the trion DOS and the effective Coulomb coupling for a 3.7 nm QD, as a function of 
the relative energy rescaled by the optical band gap. The MEG rates are fitted using averaged data 
in the manner described in Figure 2. The trion DOS is obtained using a Gaussian broadening of 10 
meV. For the largest dot (4.3 nm), the MEG rate and effective Coulomb coupling are only shown 
for electrons due to the significant computational cost required to calculate MEG rates from the 
large set of possible positive trion configurations. 
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which exhibit large values of \W\ 2 even when their wave functions are largely localized near the 
surface of the QD (Fig. 4b). Fig. 4 suggests that there is no obvious correlation between \W\ 2 
and the core/surface projection of the wave functions. Surface localization of the wave functions 
does not appear to have a dramatic effect on the MEG rate. The dependence of the MEG quantum 
yield on surface passivation, reported e.g. by Beard et al., 18 appears to be related to the effects of 
surface passivation on competing decay channels for the photo-excited electron-hole pair. 
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Figure 4: (a) Surface projection fraction, S,, of the single -particle wave function norm plotted as 
a function of the electron and hole energy. Each point corresponds to a single-particle eigenstate. 
The surface region is defined as the set of points closer than R\y$ (the Wigner-Seitz radius of bulk 
CdSe) from at least one surface Cd/Se atoms. The effective Coulomb coupling, \W\ 2 , is shown 
(solid line) for the region where MEG occurs, (b) Charge density distribution (arbitrary units) of 
hole states between -8.0 eV and -7.5 eV of a 3.7 nm CdSe QD, plotted on a (001) plane. 



Having discussed the MEG rates of individual electrons and holes, we now consider the MEG 
rate of a single exciton decaying into a biexciton. Because many-body selection rules forbid the 
simultaneous transition of the electron and the hole (in first-order perturbation theory), the MEG 
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rate of an electron-hole pair (h, e) is simply the sum of the MEG rates of the electron and the 
hole: F(h, e) = + T e , where we have assumed that (T e ) is weakly perturbed by the presence 
of a spectator electron (hole). When the QD is excited by a laser pulse, different e-h pairs that 
are nearly degenerate in energy can be excited. Thus, we calculate the MEG rate for a photon of 
energy, ha, as 



r(fto) = £p e (ftffl)r e +£/> A (ftffl)r A (i) 

e h 

where 

phi*®) = r^r'.^r (2) 



Lc\M h , c 


2 8(E h , c 


-ha) 


Lv,c\M v ,c 


\ 2 8(E V , C 


-ha) 


Iv \Mv,e 


2 S(£ v , e - 


-ha) 


Lv,c\M v ,c 


\ 2 S(E v ,c 


-ha) 



Pe(ha) = ^ J~ s r: Z\ (3) 



Here p% (p e ) is the probability of creating a hole in the valence state v (an electron in the conduction 
state c) as a result of the absorption of a photon of energy ha. The integral M V)C = (V / v| | *|Vc) is the 
dipole matrix element that characterizes the transition from the valence state v, to the conduction 
state c, and E V)C = £ c - £ v - J VjC . 

Figure 5a shows the MEG rate for all four CdSe QDs as a function of the photon energy, ha, 



calculated using [equation] [1][]1 As expected (see e.g. Fig. 2), the MEG rate increases as the 
photon energy increases. For a given absolute photon energy, the MEG rate increases as the dot 
size increases (Fig. 5a). However, the trend is just the opposite if the photon energy is rescaled 
by the optical band gap of each QD (Fig. 5b). This is the result of quantum confinement which 
causes the optical gap to increase with decreasing QD size, which implies that the photon energy 
has to increase to keep the ratio ha/E gap constant. The optical band gap of the four CdSe QDs 
in this work decreases from 2.14 eV (2.9 nm diameter) down to 1.86 eV (4.3 nm diameter). The 
calculated MEG rate is relatively low, e.g. T < 10 -2 ps _1 at ha = 2JE gap for the 3.2 nm QD. 
This result suggests that, at low photon energies, MEG can be easily overcome by competing 
carrier relaxation processes, in agreement with the experimental results of Nair et alJ^ 



10 




Figure 5: MEG rate [ [equation] [1][]1 1 as a function of (a) absolute photon energy and (b) relative 
photon energy scaled by the optical band gaps for four CdSe QDs. Dashed lines are calculated 
including only the contribution from excited electrons, i.e., Y e in | [equation] [1][] 1] For the largest 



dot (4.3 nm), the MEG rate includes only the electron contribution, as the hole contribution is 
negligible (see text). 



Figure 5 also makes clear that the main contribution to the total MEG rate originates from 
electrons decaying into negative trions. The reason is that, because of the distribution of dipole 
matrix elements [[equation] [2] [] 2 and [equation] [3][]3[ , the excess energy of the photo-generated 
electron is larger than that of the hole, so the electron has much higher probability of decaying by 
MEG, despite the fact that the MEG rate for an individual hole is larger than that of an electron 
with the same excess energy (Fig. 2). If an excited electron could be efficiently converted to an 
excited hole through an Auger process, though, the total MEG rate could be significantly enhanced. 

The trends in size scaling can be made more precise by fitting the results to power law functions 
of the QD volume. These fits are given in Fig. 6. We find that, for a given absolute photon energy, 
the MEG rate scales as Y = Tq + aV^ with /3 ranging from -0.70 to -0.91 in the energy range 
under consideration. In contrast, the MEG rate, as a function of the relative photon energy, scales 
as r = aVP with /3 lying between -1 .36 and -1.21 up to 2JE gap (Table 3). It has been argued^that 
the MEG rate should not exhibit a significant volume dependence because of the offsetting volume 
dependencies of Coulomb coupling and trion density of states. Specifically, it was suggested that 
the average Coulomb coupling scales as V~ 3 while the trion DOS scales as V 3 . Our calculations 
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Figure 6: Volume dependence of the MEG rate of CdSe QDs, calculated (a) at fixed absolute 
photon energy and (b) at fixed relative photon energy scaled by the optical band gap. Solid lines are 
power law fits to the MEG rate for prescribed energies, and the resulting coefficients are collected 
in Table 3. 



show, however, that the volume dependence of the Coulomb coupling is roughly V~ for excited 
electrons but closer to V ~ 4 for holes. The trion DOS scaling depends on photon energy and is V a 
with a 3.1 to 4.1, while for the excited holes a ranges from 3.8 to 5.5. These results show that the 
volume scaling contributions to MEG rate from the effective Coulomb coupling and trion DOS do 
not cancel each other and there is an energy dependence that must be considered as well. The net 
scaling of the MEG rate is quantified in Fig. 6 and Table 3. 

Table 3: Fitted parameters for the volume dependence of MEG rate: (left) r = Tq + aV^ at given 
absolute photon energy, h(o; (right) T = aV^ at given photon energy scaled by the optical band 
gap of the QD, h(o/E gap , as shown in Figure 6. 



heo(eV) /3 


h(o/E gap /3 


4.5 -0.70 
4.8 -0.78 
5.0 -0.82 
5.2 -0.91 


2.2 -1.36 
2.4 -1.21 

2.6 -1.20 

2.7 -1.21 



While the choice of plotting the MEG rate versus absolute or relative photon energy (see e.g. 
Figs. 5a and 5b) may appear arbitrary, it has been argued that for a given QY the MEG efficiency 
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increases as the excitonic gap increases, because the energy of the multi-excitons produced by 
MEG increasesJ 2 -^ 2 ^ Beard et al.^have suggested that the MEG efficiency is proportional to the 
slope of the QY vs. the relative photon energy h(o/E gap (above the MEG threshold h(O t h). Thus, 
the dependence of the MEG rate on QD volume calculated at fixed relative photon energy (Fig. 
6b) is more appropriate to describe the trend of the MEG efficiency with QD size than the MEG 
rate calculated at absolute photon energy (Fig. 6a). 

In summary, we have applied an atomistic pseudopotential method to quantify the relation 
between size and MEG rate in CdSe QDs. The MEG rate increases as QD size decreases only if 
the comparison is made for the same relative photon energies, i.e., after normalization using the 
volume dependent optical gap. The opposite trend is found, however, if the comparison is made 
using the same absolute photon energy for each QD. Both the effective Coulomb coupling and the 
trion DOS are energy dependent with the latter being a stronger function of energy. Although the 
decay rate of an excited hole is much larger than that of an excited electron, the total MEG rate is 
dominated by photo-excited electrons, because they have higher excess energy than photo-excited 
holes. Volume scaling laws for overall MEG, trion DOS, and Coulomb coupling integral have now 
been quantified under realistic conditions and with a minimal number of simplifying idealizations. 

Methods 

Our computational approach utilizes an atomistic pseudopotential that is tailored to capture the 
experimental electronic structure of CdSe QDs.^^ The CdSe quantum dots are constructed by 
carving out a Se-centered sphere from bulk CdSe zinc -blende structure. The bulk crystal utilizes 
the experimental lattice constant of ao=6.08 A.^ All surface Cd and Se atoms are passivated with 
ligand-like atomic pseudopotentials that are fitted to remove the surface states from the band gap 
of the QD. The Cd and Se atomic pseudopotentials are fitted to experimental bulk CdSe properties 
such as transition energies and effective masses and to single-particle bulk wavefunctions calcu- 
lated from first principles J 28 * 29 ! The single-particle energies, £,-, and wave functions, i/a,-, are then 
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obtained by solving the single-particle Schrodinger equation 

h 



l-^V 2 + V(r)+V S0 ]V l (r,o)=£ l \ l r l (r 1 o) (4) 
2m 



Here the single particle wave functions, i/a,-, are expanded in a plane wave basis set, Vso is tne 
nonlocal spin-orbit operator, and the local potential, V(r), is a linear superposition of screened 
pseudopotentials centered at the atomic positions, {if}, of atomic type a: 

V(r) = J>a(r-*) (5) 

The impact ionization rates for electrons and holes are calculated here by applying Fermi's golden 
rule to the formation of trion states: 

r e = ^\(e\W\f)\ 2 8(Ej-E e ) (6) 
n f 

r h = -Z,\(h\W\f)\ 2 8(E}-E h ) (7) 
n f 

Here \e) and \h) are the initial, excited states for an electron in the conduction band and a hole in 
the valence band, respectively. The operator, W, is the screened Coulomb matrix element between 
the initial and final states. The final trion states are |/) = |«i,C2,^i) and |/) = \e\,h\,h2), 
respectively, and their energies are calculated as 

Ef = (£ e l +£e2~ + (Jel,e2 — h\,h\ ~ hl,h\)i (8) 



and 



E~c = (E e i — Ehl — Efa) + (J hi, hi — J el, hi —JelJo)- (9) 
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The screened Coulomb interactions, J t -j are given by 



e r.r r — r 



7TlV0'( r ' CT ) 



(10) 



where e(r, r') is the dielectric screening function of the QD.^We adopt a Lorentzian line shape, 

in evaluating 



[equation] [6][]6 



and 



[equation] [7] []7 



to account for 



phonon broadening effects. A value of 7 = 10 meV is used for the summation over the final 
states; 34 changes in 7 from 5 meV to 20 meV were found to not meaningfully affect MEG rates. 
Rabani et al.^ used a pseudopotential approach to calculate the MEG rates of small (up to 500 
atoms) CdSe and InAs nanocrystals. Our approach builds on their method by (i) including spin- 
orbit coupling, which is significant in CdSe quantum dots, and (ii) including excitonic interactions 
[equation] [6] [] 6| and |[equation] [7] []7 Electron-hole Coulomb interactions |[equation] [ 10] [] 1 0| 



m 



lower the energy of the final (trion) states [see [equation] [8] [] 8 and [equation] [9] []9|, thereby 
increasing the density of final states which are in resonance with the initial state. The net result is a 
significant increase in the MEG rate (an order of magnitude) compared to the case where Coulomb 
interactions are neglected. 
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